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Abstract 

We study the influence of boundary conditions on stationary, periodic patterns in one-dimensional 
systems. We show how a conceptual understanding of the structure of equilibria in large domains can 
be based on the characterization of boundary layers through displacement-strain curves. Most promi¬ 
nently, we distinguish wavenumber-selecting and phase-selecting boundary conditions and show how 
they impact the set of equilibria as the domain size tends to infinity. We illustrate the abstract 
concepts in the phase-diffusion and the Ginzburg-Landau approximation. We also show how to 
compute displacement-strain curves in more general systems such as the Swift-Hohenberg equation 
using continuation methods. 
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1 Introduction 

Stripe patterns are arguably the simplest non-trivial patterns observed in nature. Stripes form in 
numerous contexts, ranging from classical fluid experiments such as Taylor vortices in the Taylor- 
Couette flow and convection rolls in Rayleigh-Benard convection, over bilayers in diblock copolymers, 
to granular media, and reaction-diffusion systems. A well-studied context of stripe formation is an 
instability of a spatially homogeneous state in a spatially extended medium. A linearized analysis can 
then predict wavenumbers as linearly fastest growing modes. While known for more than a century in 
fluid dynamics, the occurrence of patterned states as fastest growing modes was conjectured by Turing 
in 1952 [18] for reaction-diffusion systems and experimentally realized in 1991 |3|. 

The arguably simplest example for the formation of stripes is the celebrated Swift-Hohenberg equation 

ut = —(A -b l)^u + ftu — u^. 

Considered on x G M, for instance, the linearization at the trivial state u = 0, 

Ut = -(A -b l)^u -b ftu, 
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Figure 1.1: Boundary layers selecting different wavenumbers and phases of the periodic pattern Ugt (top left); 
boundary layers for left and right boundary, ui/r, matched at x = 0 (bottom left). Schematic plot of boundary 
manifold B, strong stable fiber IF®® intersecting B and thus yielding a boundary layer, the family of periodic 
orbits, and the unstable fibers (right). 


can be readily analyzed using Fourier transform, 

^u{k) = (-(1 - + n)u{k). 

For /i > 0, wavenumbers with \k‘^ — 1| < are unstable. The fastest-growing wavenumber is |/c| = 1. 
In fact, for all fi < 1, only wavenumbers with k ^ 0 are unstable, so that for small amplitude only 
spatially patterned perturbations are amplified. 

A nonlinear analysis readily reveals that together with the linear instability, nonlinear patterns bifurcate 
from the trivial state. For /j. > 0 small, one finds 27r//c-periodic “striped” solutions Ust{kx; k, fi), 
Ust{^;k,^) = Ust{^ + 27r;A:,/i), for all unstable A:-vectors, \k‘^ — 1| < Many of those periodic 

solutions are compatible with the boundary conditions in large but finite systems. With, say, periodic 
boundary conditions and system size L, we need to require k G (27r/L)Z, which yields O(Ly^) different 
A:-values. 

For boundary conditions other than Neumann, Dirichlet, or periodic, equilibria are difficult to char¬ 
acterize completely. In fact, bifurcations diagrams tend to be very complex, and of little help when 
trying to “describe” typical dynamics. Interestingly, different types of boundary conditions can provoke 
quite different types of dynamics. For instance, heated walls in Rayleigh-Benard convection are known 
to favor parallel alignment of roll solutions with the boundary, as opposed to the otherwise typical 
perpendicular alignments. Slow drift of rolls has also been observed in contexts of Rayleigh-Benard 
convection. 

The purpose of this work is to study the role of boundaries in such pattern forming systems in a 
systematic fashion. The key idea is to decompose the problem of describing equilibria in large finite 
domains into two steps (see Figure o for an illustration): 

(i) describe equilibria in semi-infinite domains x G M+ and x G M-; 

(ii) patch equilibria from semi-inhnite domains at x = 0 to obtain equilibria in bounded domains 
X G (—T, L'j . 

The boundary conditions in step (ii) at x = L should be studied as boundary conditions at x = 0 in 
step (i) for x G M_, boundary conditions at x = —L give boundary conditions for x G M+. 
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Outline. We study a general class of equations in Section In particular, we characterize stripes, 
introduce spatial dynamics, define boundary layers, and discuss displacement-strain relations for bound¬ 
ary layers. We also present our main general matching result for bounded domains. Section [^illustrates 
the concepts at the hands of the (somewhat trivial) example of the phase-diffusion equation. In Sec¬ 
tion [^ we study the specific (integrable) example of the Ginzburg-Landau equation, which arises as 
a modulation equation for small-amplitude striped patterns. Section outlines an effective numerical 
procedure for the computation of displacement-strain relations and illustrates results for the Swift- 
Hohenberg equation. We conclude with a discussion in Section 


2 Boundary layers for striped phases: wavenumber and phase selec¬ 
tion 


We define families of stable stripes in Section \2.1\ and introduce spatial dynamics in Section |2.2[ We 
introduce boundary layers in Section |2.3[ We briefly explore stability in Section 2.4 Section 2.5 


IS 


concerned with gluing boundary layers to obtain equilibrium configurations in large but finite domains. 
Section 2.6 formalizes the gluing process as a subtraction of curves followed by quantization. Section 


2.7 contains a list of conceptual examples. 


2.1 Families of stripes 

We consider a general semilinear parabolic system on the real line 

ut =+ f{u,dxU,... u£R^,x£R (2.1) 

where / is smooth and reflection symmetric, 

/(u, Ml, U2 ) • • • ) U2m—l) — /(^j ) '^2i ■ ■ ■ i '^2m—1) • 

Such systems may possess even, periodic solutions Ust, which we refer to as stripes. The linearization 
at such stripes, 

2m 

Cstv = -{idxf'^v + ^ djF{u, dxUst, ■■■, 
j=i 

defines an elliptic operator on L‘^{R,R^). The spectrum of L is the union of the point spectra of the 
associated Bloch operators, 

2m 

= -{i{kd^ + 'd)f'^v + '^djf{u,dxUst,. ■. ,dl'^~^Ust){kd^ + 'Fy~^v, 

i=i 

where the operators are defined on L^((0, 27r), M'^) with periodic boundary conditions. Since 
u'g^{kx]k) contributes to the kernel of Cq, one always finds a band of eigenvalues X{i), A(0) = 0, to 
Ci- Symmetry shows that A(^) = G M. We therefore say that a periodic pattern is stable if 

the spectrum of C is contained in Re A < 0, except for the critical branch of eigenvalues A(£), for 
which we require A''(0) < 0. Such spectral stability has been shown to imply nonlinear stability with 
respect to localized perturbations and we will henceforth refer to this criterion simply as “stability”; 
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see [HlElIle]. Instabilities can occnr due to other bands of spectrum crossing the imaginary axis, or 
due to an Eckhaus instability, when A"(0) changes sign [TT] , 

Hypothesis 2.1 (Families of stripe solutions) We assume that there exists a bounded family of 
stable stripes Ust{kx; k),for wavenumbers k ^ Jk = (^min, ^max); 0 < fcmin < ^max < oo, 

Ustif,-, k) = Ustii + 27r; k), Ustit k) = k) Ust{-] k) G C‘^{Jk, <^^^(0, 27r)). 

We emphasize that the interval Jk is not assumed to be small. In the examples that we consider below, 
the boundaries A^min/max determined by Eckhaus instabilities. Most of our analysis is insensitive to 
oscillatory instabilities, A G iM \ {0} so that in this case the effective range of wavenumbers could be 
further extended. 

Lemma 2.2 (Robustness) Given a stable stripe, there exists a family of stripes nearby. Moreover, 
families of stripes depend smoothly on system parameters. 

Proof. The kernel of the Bloch-wave operator Cq is trivial in so that one can continue a 

stable periodic pattern smoothly in the period. Since the condition A''(0) < 0 is open, all members in 
the local family will be stable. One similarly establishes persistence of global families with respect to 
parameters. ■ 

Remark 2.3 (Minimal period) We may assume throughout that the period 2 tt is the minimal period. 
We may choose the period to be minimal for a fixed member in the family. Changes of the minimal 
period always imply a bifurcation, hence would imply additional spectrum at X = 0 and contradict our 
assumption on stability. As a consequence, the stripes have precisely two reflection symmetries, x —>• —x 
and X —)• 27r — X. 

2.2 Spatial dynamics 

The steady-state equation 

0 = + f{u, d^u,..., 

can be rewritten as a first-order system, 

= F([/), U = {u,d,,u,...,dl'^-\) (2.2) 

which possesses a reversibility symmetry 

F{RU) = -RF{U), R{uo, Ul, . . . , U2m-l) = (m, -Ul, . . . , -U2m-l)- 

Stripes correspond to reversible periodic solutions RU{kx;k) = U{—kx\k). Purely imaginary Floquet 
exponents of the linearization at such a periodic solution give rise to kernels of the operators in 
a straight-forward fashion. Moreover, the algebraic multiplicity of the simple zero-exponent is given 
by the order order of tangency of the curve of critical spectrum A(^); see for instance [141 Lemma 
2.1]. The family of stripes therefore forms a two-dimensional normally hyperbolic manifold S in 
Reversibility implies that both stable and unstable manifold are mN + 1-dimensional, and 

mapped into each other by R. Both are smoothly fibered by mN — 1-dimensional strong (un)stable 
fibers ly/"". 
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2.3 Boundary layers 


We consider (2.1) on x G M+, supplemented with mN boundary conditions, which we write in the form 
[/(x = 0) G where B is an mWdimensional smooth manifold in Natural examples are 


• Dirichlet, di^u = aj,0 ^ j < mN] 

• Neumann, dx^~^^u = aj,0 ^ j < mN; 

• clamped, diu = aj,0 ^ j < mN; 

• free, diu = aj,mN ^ j < 2mN. 


The terms clamped and free allude to the plate equation. Of course, many mixed and nonlinear 
variations are possible, and not all guarantee well-posedness of the parabolic equation ( 2 . 1 ). 


Equilibria of (2.1) in 




that satisfy the boundary conditions B are bounded solutions 


U{x) to the ODE (2.2), with U{0) G B, where B is now considered as a natural mAi-dimensional 
submanifold of the phase space of of (2.2). In other words, equilibria are naturally identified 

via a multiple shooting problem; see Eigure [m Equilibria of interest to us converge to stripes as 
X —)> oo. They therefore lie in the intersection B 0 1T“. Simple dimension counting suggests that this 
intersection is one-dimensional. Thom transversality indeed guarantees just this for generic manifolds 
B. Moreover, the intersection projects smoothly along the fibers onto a curve in S. Parameterizing S 
naturally by /c G and an angle ip, we can now characterize boundary layers through curves in the 
cylinder ip, k £ x J^. 


Definition 2.4 (Boundary layers and displacement-strain relations) We say a boundary layer 
UhiiO *-5 transverse ifTjj^^p;p^B fh is one-dimensional, not contained in . We refer 

to the base point of the fiber IT®® that contains ^bi(O) as the asymptotic phase ip and wavenumber k. 
If all boundary layers are transverse, we refer to the collection of base points as displacement-strain 
{k, ip) curves 7 . Unless otherwise noted, we parameterize displacement-strain curves by arc length j(s). 
ITe refer to a relation d{k, ip) which vanishes precisely at 7 as displacement-strain relation; see ^2^ 
for an explicit formula. 


Remark 2.5 • The terminology displacement-strain refers to the fact that the asymptotic stripes 

are “displaced” by ip relative to a fixed, even, stripe pattern on x G M. Strain refers to the intuitive 
compression and expansion of the stripe relative to a fixed wavenumber. 

• The definition of ip in S is unique with the convention S = {Ust{kx — ip;k), k G Jk}, when 
requiring Ust{-,k) to be even and after fixing f7st(s^o) (which rules out a shift by n). We then 
have 

\U{x)-U,t{kx-ip;k)\^Ce-^^, (2.3) 

for some positive constants C, rj. 

• In a completely analogous fashion, one can define boundary layers for right-bounded domains, 
X G M-, 17(0) G B-\., and 

\U{x)-U,t{kx-ip;k)\^Ce^U (2.4) 
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Note that boundary layers in left-bounded domains give boundary layers in right-bounded domains 
via reflection, but for reflected boundary conditions RB+ = B-, replacing dx by —dx- 

• We emphasize that the curve 7 may have self-intersections, stemming from the projection of a 
smooth non-intersecting curve along the strong stable fibration. 


Given the (generic) transversality in Definition 2.4 we can envision families of boundary layers, globally 
in k, Lp. Note however that even when all boundary layers lie in a bounded set in phase space , the 
resulting displacement-strain curves may still terminate, for instance when the p — fc-curve approaches 
the boundary of Jfc,or when the “length” of the heteroclinic diverges; see [ 6 ] for some context on 
homoclinic and heteroclinic continuation. We will in fact encounter both possibilities in the example 
of the Ginzburg-Landau equation. Section 


Definition 2.6 (Phase- and wavenumber selection) A displacement-strain curve 'y C x 
encodes 


• wavenumber selection, if Pk'~i / Jki 

• phase selection, if . 

Here, P^/k denote projections onto the ip- or k-component, respectively. If ^ is a closed curve, we write 
i{'y) for the winding number of the curve P^y in . 


We will commonly encounter winding numbers 0 and 1, but also an example if winding number 2 in 
the Ginzburg-Landau example. 

We will explore a definition of boundary layers not based on spatial dynamics in Section This 
definition will prove useful when constructing boundary layers numerically. 


2.4 Stability of boundary layers 

Without any further assumptions, little can be said about the stability of boundary layers. Given the 
stability of the asymptotic stripe pattern, we can say that the essential spectrum is contained in the 
negative half plane with the exception of a branch of essential spectrum touching the imaginary axis; 
see [3]. The mere existence does, of course, not give information on possible oscillatory instabilities, 
but we do have information on possible zero eigenvalues. We therefore count the parity of unstable 
eigenvalues of the linearization at a boundary layer Pbi = 

2m 

Chiv = -{idxf'^v + ^ djF{u, dxuu, ■ • •, dl'^~^uu)di~^v, x > 0, 

equipped with boundary conditions {v, dxV, ..., G TB. In other words, p = 1 refers to an even 

(or zero) number of unstable eigenvalues of £bi and p = — 1 to an odd number. 

The following lemma relates the infinitesimal displacement-strain relations to stability. 
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Lemma 2.7 (Parity and displacement-strain) For any family of boundary layers and correspond¬ 
ing displacement-strain curve 7(5), there exists a parity e G {±1} so that 

p = e - sign Pk 

Proof. We need to show that A = 0 belongs to the extended point spectrum precisely when k' = 0 
along a family of boundary layers, and that we have a strict crossing whenever k' changes sign. 

Under our assumption on symmetry and stability of periodic patterns, the essential spectrum in the 
neighborhood of the origin is given by the negative real axis A ^ 0. We consider the linearized equation 
at a boundary layer, with spectral parameter A = 7 ^. Since coefficients converge exponentially, there is a 
complex linear change of variables, analytic in 7 and smoothly depending on x, so that the strong stable 
{Nm — 1 -dimensional), strong unstable (Nm — 1 -dimensional), and center subspaces 2 -dimensional are 
independent of x and 7. in these coordinates, the (linearized) boundary conditions at x = 0 define 
an A^m-dimensional subspace that depends analytically on 7. We list vectors in these coordinates 
as ([/^, U®®, [/'^'^)^. Since 7 = 0 is a simple branch point of the dispersion relation, we can choose 
coordinates in the center subspace such that the stable subspace is spanned by the Nm — 1 column 
vectors in (0,id, 0)^ and the vector ((1, —7), 0,0)^ -|- 0(7^). The boundary conditions are spanned by 
the Nm—1 column vectors in (0,0, id)^ and (Q^, 0,0)^, where Q is an invertible linear transformation. 
We can form the determinant £{'j) from basis vectors of stable eigenspace and boundary tangent space 
to build the Evans function which then tracks eigenvalues near the origin 7 = 0 . Eigenvalues in 7 > 0 
correspond to actual unstable eigenvalues, whereas 7 < 0 corresponds to resonance poles. Since the 
stable subspace limits on the eigenspace at 7 = 0 , which is spanned by the derivative of the periodic 
pattern, we find that Q'j' = {qu(p' ^22^0^i with qii,q 22 7 ^ 0 from invertibility. We therefore 

find that 

£ (7) = (922^' - qiiF'i + 0(7^)) ^0(7), 

where £/ is analytic and nonzero near the origin. We conclude that 7 = 0 is a root precisely when 
k' = t), and the sign of 7 changes precisely when the sign of k' changes. ■ 



Remark 2.8 (Nonlinear stability) Given the results on stability of stripes on M, one expects spec¬ 
trally stable boundary layers to he stable. In fact, one would expect a slightly faster pointwise decay, 
~ given that the Evans function does not vanish at the origin for boundary layers, but does vanish 

for stripes. This is in analogy to the heat equation with, say, Dirichlet boundary conditions on the half 
line. 




Stable 


Remark 2.9 (Orientation) We can orient displacement curves so that p = sign 
boundary layers within a connected family then necessarily either compress or expand along 7 . The 
somewhat trivial examples of phase diffusion, below, illustrate that there is no obvious way to avoid an 
ambiguity in orientation (or, equivalently, choice of e in Lemma 2.1). 


Remark 2.10 (Displacement-strain dependence) The results on stability suggest that, typically, 
k could be thought of as a parameter so that turning points correspond to exchange of stability as in a 
traditional saddle-node bifurcation. It is mainly this motivation that leads us to displaying displacement- 
strain relations with k as the independent, horizontal variable in the remainder of this paper. 
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2.5 Finite domains — gluing boundary layers 


In this section, we describe equilibria in large but finite domains, based on left and right boundary 
layers. The key conclusion is that displacement strain curves for boundary layers associated with 
boundary conditions on the left and right side of the domain allow for a complete description of the 
set of equilibria for large enough domains {—L,L), with boundary conditions at ±L. . 

We employ a two-sided shooting approach. We start with transverse boundary layers Ui{x] k, ip),x > 0 
and Ut:{x; k,(p),x < 0, which satisfy the boundary conditions at Bi/j. at x = 0, and 

Ui{x] k, (p) — Ust{kx — (p; k) —)-0, x —)• -|-oo, 

Ur{x] k, (p) — Ust{kx — (p; k) —)-0, x —)• —oo, 

with associated displacement-strain relations 

d\/rik,ip) = 0 . 

Shifting C/pj. to the boundaries, we look for solutions U that are close to U\{x + L-,k\,(p\) on (—L,0) 
and Ur{x — L;kr,p>r) on (0,L). Continuity of the solution at x = 0 up to exponentially small terms 
gives phase and wavenumber matching, 

kiL — ipi = —krL — (pr mod 27r, k\ = kr- (2-5) 

The following proposition shows that this formal phase and wavenumber matching procedure actually 
determines solutions. 

Proposition 2.11 Assume the existence of families of transverse boundary layers U\{x-,k,ip),x > 
0 and Ur{x;k,(p),x < 0, to boundary conditions Bxj^ at x = 0, with displacement-strain relations 
d\i^{k,(p) = 0. Then there exist solutions on {—L,L) with boundary conditions Bxj^ at x = ±L which 
are L°°-close to Ur{x — L; k\, ip\) on x < 0 and to U\{x -\- L; k^, p>r) on x > 0 if exponentially corrected 
phase and wavenumber matching conditions hold, 

k\ = kr TZk (2-6) 

{ki -\- kr)L = ipi — ipr 7^,^mod27r, (2-7) 

where TZj{k\, kj-, tpi, ‘Pr, L), j = p, k, are smooth and TZj = 0(e“^^), together with its derivatives. 

Proof. We write U = Ui{x + L)y_(x) -|- Ur{x — L)x+{x) + W{x), where x± is a smooth partition 
of unity, x_ -|- x+ = 1, x±{x) = 1 for ±x > 1. Exploiting the fact that the boundary layers solve the 
ODE, we find the equation 

W = F(Ui(x + L)x-(x) + Ur(x-L)x+(x) + W(x))-F(Ui(x + L)x-(x) + [/r(x-L)x+(x)) + R(x), (2.8) 
with boundary conditions 

W{±L) G Bx/r - t/i/r(0), (2.9) 

where the remainder is 

R{x) =F{Ux{x + T)x-(x) + Ur{x - L)x+(x)) - F{Ux{x + T))x-(x) - F{Ur{x - L))x+{x) 

- Ux{x + L)x'_(x) - Ur{x - L)x+(x). 



We view the system (2.8) as an equation with variables W and r±, where r± enter as a parameterization 
of the boundary layers 


Ul/r(x) = Ui/j.{x; ki/^{T±), </3l/r(r±)), dl/r(fcl/r(r±), VPi/r(rV9i/r(r+))) = 0, 


Note that the remainder is small when the system (2.6)-(2.7) is satisfied, 
stand the linearized operator 


We therefore need to under- 


W - F'([/i(x + L)x-(x) + [/r(x - L)x+(x))W, W{±L) G TBy,\uyAoy 

Note that this operator is close to the linearization at individual boundary layers on {—L, 0) and (0, L), 
respectively. This linearization at individual boundary layers is not invertible in the limit L = oo. 
It is however Fredholm in exponentially weighted spaces, with index -1. We therefore study relaxed 
boundary conditions 

W{±L) + ai/r'0i/i.(O) G Bi/j. - Ui/riO), (2.10) 

. Here, aj/r are independent scalar variables and ^jJl/^{x) are solutions to 

V’l/r - = 0, 

such that '01/1. 0 THi/r, with at most linear growth at ±oo, respectively. 

We will show that the associated linearized operator, 

W - F'([/i(x + L)x-(x) + Ur(x - L)x+(x))W, W(±L) + ai/r0i/r(O) G 

is invertible with L-uniform bounds in exponentially weighted spaces, 

\\W\\l2 = ||min{e^("=+'^),e’'(-^+^)}iy(x)||i2, (2.11) 

with domain where IF, W' G L^. 

For this, we solve = H hy decomposing H = + i7+, H± = x±H, 

W. - F'(Ui)W. - [F'(UiX- + UrX+) - F'(Ui)jW+ = H_, 

W+ - F'(Ur)W+ - [F'iUiX- + UrX+) - F'(Ur)jW. = H+, 

1F_(-L) + lF+(-L) + a_V’i G TBi, 

W+{L) + W-{L) + a+0r G TBr. 


replacing (2.9 


In the exponentially weighted norms (2.11), the boundary terms W+(—L) and W-{L), 
coupling terms involving W+ in the first equation and W- in the second equation are 
small, so that this system is a small perturbation of 


as well as the 
exponentially 


W- - F'{Ui)W- = H_, 
W+ - F'(Ur)W+ = H+, 
1T_(-L) + a_0i G TBi, 
lF+(L) + a+0r em. 


We can extend F[± to half lines and thereby find W± by inverting the linearization at the boundary 
layers. 
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Summarizing, we have shown uniform bounded invertibility of the linearization of (2.8), with extended 
boundary conditions (2.10), variables W, a±, and parameters r±, L. The residual is small, R = 0(e~^^) 
when the phase matching conditions. 


ki(T_) = kr(T+), 2krL = V?l(r_) - (Pr{T+), 

are satisfied. We therefore obtain families of solutions (W,a±) = {W*, a^){T±, L). It remains to solve 

a+(r+,r_,L) = 0, al(r+, r_, L) = 0. 

A straightforward but tedious expansion of this system shows that it is equivalent to wavenumber and 
phase matching, up to exponentially small corrections, 

k] = kr + k\L — Lp\ = k^L — (pr + mod 27r, 


where 

{kup\) = iki{T-),(pi{T-))+ a-{k[{T-),p[{T-))-^. 
This establishes the desired phase and wavenumber matching conditions. 


Remark 2.12 A more geometric proof would he based on flowing the boundary conditions from x = ±L 
to X = 0 and exploiting inclination lemmas (see for instance m Thm 3.1]) to see that the interseetion 
is transverse in a complement of the tangent space to the two-dimensional family of periodic orbits, so 
that intersections can he computed from matching within this manifold up to exponentially small terms. 
Our approach here is closer in spirit to m- 


2.6 Subtraction and quantization — from displacement-strain curves to equilibria 

In the following, we discuss conceptually the solutions to (2.6)~(2.7). We will neglect exponentially 


small terms, which can be easily accommodated in the analysis. As in the proof of Proposition 2.11 
we parameterize solutions to the displacement-strain relation as curves via 

rpr = {{h/rij±),ip\/flT±), 


with nonvanishing tangent vectors \{'p',k')'\^ = 1 and study wavenumber and phase-matching, (2.6) 


(2.7) as two equations in three variables t±,L. It is conceptually helpful to add an intermediate step. 


solving 


0 = fci(r_) — k, 

0 = kr{T+) — k, 

0 = (^i(r_) - Pr{r+) - Ip- 

Note that with the definition of k,'if from this system, phase matching reduces to ip = 2/i;Lmod27r. 

The linearization of this auxiliary system is readily seen to be onto as long as k[ and fc' do not vanish 
simultaneously. Since in this case the linearization with respect to k, p, and either r+ or r_ is invertible, 
we can parameterize the solution set as curves (fc, V’)('r). Critical points of A:(r) are given by the union 
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Figure 2.1: The figures show schematic plots of left and right displacement-strain relations, their differential 
displacement-strain curves, and the quantized intersections with ip = 2kL. From left to right: two-sided phase 
selection, phase and wavenumber selection, two-sided wavenumber selection with Ji C Jr (note the two branches 
of the differential displacement-strain curve), and two-sided wavenumber selection with marginal overlap. 



Figure 2.2: The figure illustrates the effect of increasing the domain size on equilibria. Wavenumber selection 
with nonzero winding number (left) leads to a continuous (snaking) family of equilibria as L is increased. Phase 
selection (right) leads to chains of creation and annihilation events of equilibria near the Eckhaus boundary. 


of critical points of ki/j.{T±). Critical points of correspond to values of k where tangent vectors of 
left and right displacement-strain relations are collinear, 

k{ K 
i’i V'r 

We can formalize these operations as follows. 


Ip' = 0 


Definition 2.13 (Differential displacement-strain) We refer to the curves {k,'il>){T) as differen¬ 
tial displacement-strain curves, obtained by taking the difference between left and right displacement- 
strain curve, pointwise in k. 


From differential displacement-strain curves, we obtain equilibria in the bounded domain for large L 
by intersecting with the straight line if = 2kLmod 2 tt, effectively quantizing the curves. 

In summary, the operations that lead to a description of equilibria in bounded domain are: 


• vertical subtraction, which gives differential displacement-strain curves; 

• quantization of differential displacement-strain, by intersecting with = 2kL. 


2.7 From displacement-strain to eqnilibria: conceptnal examples 


In the following, we analyze several representative cases in some more detail; see Figure 2.1 for illus¬ 
trations of differential displacement-strain curves and Figure [2]^ for the effect of varying L in the case 
of wavenumber and phase selection by the differential displacement-strain curves, respectively. 
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Two-sided phase selection. In this case, both displacement-strain curves are defined for k^-ji G J, 
and so is the differential displacement-strain curve. In the simplest case, if) = V’(^) is a graph, and 
intersections with the straight lines 2kL are monotone in L, increasing from k- to A:+ as L increases. In 
particular, specific equilibria appear and disappear as L is increased by an order-one amount because 
of instabilities at the boundary of J. 

Phase and wavenumber selection. Suppose that k^ G J and ki G Ji. The differential displacement- 
strain curve is now defined on Jy, as a closed curve, with the same winding number of the V' component 
as (^r. Intersections with the straight lines 2kL give continuous families of equilibria in the case of 
nonzero winding number (the right boundary layer is only selecting wavenumbers, allowing for arbitrary 
phases). 

Two-sided wavenumber selection. When C J\, the situation is equivalent to left-sided wavenum¬ 
ber selection. New phenomena appear when and Jy are not contained one in the other. Of course, no 
equilibria exist when these intervals are mutually disjoint. In the simplest case, a mutual intersection 
gives figure-eight type differential displacement-strain curves. Again, one expects phase slips as the 
size of the domain is gradually increased. 

Summary. Envisioning slowly increasing the domain size L, we expect continuous dependence of the 
solution on L only when one of the boundary conditions is wavenumber selecting with nonzero winding 
number in the phase, and when all wavenumbers of this boundary condition are compatible with the 
other boundary condition. Other combinations will typically lead to phase slips. 

Incompatible boundary conditions (selecting disjoint wavenumber intervals) will lead to phase drift as 
we shall see in the simple example of the phase-diffusion problem, below. 


3 Boundary layers in phase-diffusion 

Using a multiple-scale expansion, one can derive simplified dynamics near roll solutions, known as 
the Cross-Newell phase diffusion equation. One expands an Ansatz Ust{0{t, x)] 6x{t, x)) assuming long- 
wavelength modulations of the phase 9, and finds the nonlinear diffusion equation 

9t = {b{9x))x, 

where b'{k) = X"{0). In particular, the long-wavelength problem is well-posed when b' > 0. Boundary 
conditions are given by smooth curves B in the (0,6*x)-plane. Each point {0,6x) = ((/^(s),/c(s)) G B, 
gives rise to a unique equilibrium in x > 0, ip + kx. Linearizing the equation and boundary conditions 
gives 

9t = b'{k)9xx, k'9 - if'9X = 0 . 

Patterns with k'ip' < 0 are unstable with eigenfunction 7 = k' j ip', and eigenvalue A = 7^. 

Stable patterns possess k'ip' > 0, with resonance pole A = 7^. Note that when ip' —)■ 0, an eigenvalue 
(or resonance pole) disappears at infinity, not unsurprisingly due to the singular perturbation in the 
boundary conditions when passing from Dirichlet to Robin (or mixed). The allowed interval J is 
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given through b'{J) > 0. Clearly, Dirichlet boundary conditions 9 = (pQ select a phase and Neumann 
boundary conditions 9x = fco select a wavenumber. In reaction-diffusion systems or the Swift-Hohenberg 
equation, phase selection is accomplished by Neumann boundary conditions. Since perturbations from 
Neumann to mixed are non-singular, we can exclude eigenvalues near infinity and changes of stability 
when (p' = Q for such more general systems. 

Within the class of phase-diffusion problems, one would require (p' to be nonzero to guarantee well- 
posedness. As a consequence, one finds k = k{p) and isolas are excluded, that is, wavenumber selection 
implies a non-zero winding number of the displacement-strain curve. 

Imposing different wavenumber selection conditions 9x = k± at x = ^L, one finds drifting patterns. 


9{t, x) 



fk- +k+ 

V 2 



dy 



t, 


which are effectively time-periodic solutions to the equation since 6 G M/27rZ. Roughly speaking, h{6x) 
interpolates the selected wavenumbers linearly, leading to a nonlinear phase profile and constant phase 
drift. Note that the expression is well-defined as long as the interval of wavenumbers defined by the 
boundary conditions lies within the stable regime, b' > 0. 

On the other hand, our analysis validates the phase-diffusion approximation in semi-bounded (or even 
bounded) domains for general striped phases in a very concise fashion: boundary conditions should be 
replaced by displacement-strain relations. 


d{e,ex) = d, 

at x = 0, respectively, where 9 is understood modulo 27r, and the equation in the interior of the domain 
should be simply 9xx = 0. The fact that displacement-strain relations may not lead to well-posed 
boundary conditions points to the fact that some defects (here, boundary conditions) can simply not 
be incorporated into a phase-modulation description. On the other hand, one can expand the effective 
boundary conditions near 9x = k^, and solutions near 9 = k^x -|- p, to obtain locally valid effective 
Robin boundary conditions. 

4 Boundary layers in the Ginzburg-Landau equation 

The Ginzburg-Landau equation 

At = Axx + A — A|Ap. A G C, 

arises as a modulation equation at the onset of a Turing instability and as such provides a universal 
model for small-amplitude Turing patterns. We will study the effect of boundary conditions within 
this equation but will also comment on the effect of scalings used to derive Ginzburg-Landau on the 
boundary conditions. Stationary solutions satisfy 

Ax = B, Bx = -A + A\A\\ (4.1) 

Exact solutions include the periodic patterns A{x) = y/l — and the defect solutions, 

Ad(x;/c) = (^y2k + i\/l — 3/c2tanh(\/l — Sfc^x/v^)^ e‘^^. (4.2) 
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One readily calculates that the argument of A{x', increases by 2arctan(-v/l — 3k‘^/{y/2k)) when 

k > 0 and decreases by the same amount when k < 0. For A; = 0, the argument is not defined at x = 0 
and the final phase shift is tt. The defects are invariant under the reflection A{x) i—>■ A{—x). Defects 
are known to be unstable for energetic reasons. 


As a Hamiltonian system with two degrees of freedom and a phase rotation symmetry, (4.1) is integrable 
and all solutions are explicit, in particular the stable and unstable manifolds of periodic orbits. 


Indeed, the bounded part of stable and unstable manifolds is given by the defect, and the unbounded 
parts are given by 

A{x;k) = (^\/2k + i\/1 — 3k‘^ coth(-\/1 — 3k'^x/\/2)J e'^*, (4-3) 

with X < 0 or X > 0, respectively. 


The Ginzburg-Landau equation arises universally as an amplitude equation for small striped patterns, 
and therefore is a prototype for our study. On the other hand, the additional symmetries in this equation 
sometimes distort phenomena. In fact, one can further “reduce” the Ginzburg-Landau equation by 
assuming slowly varying phase near periodic patterns, and then derive the phase-diffusion equation, 
discussed in the preceding section. Our main findings here do replicate the phenomena for the phase- 
diffusion equation. However, we find several new phenomena, in particular: 


• displacement-strain curves with winding number 2 (inhomogeneous Neumann); 

• displacement-strain curves where ip is not monotone (inhomogeneous Dirichlet); 

• defect pinch-off at terminal points of displacement-strain curves (inhomogeneous Neumann and 
Dirichlet); 

• reconnection crises near the Eckhaus boundaries (inhomogeneous Neumann and Dirichlet). 


4.1 Symmetries and invariant coordinates 


While one can find boundary layers using the explicit representation of stable manifolds in (4.2) and 
(4.3), we found it significantly more manageable (albeit, presumably, equivalent) to use different coor¬ 
dinates and reduce to polynomial equations. 

The stationary Ginzburg-Landau equation ( |4.1[ ) is invariant under the gauge and reversal symmetries 
n^{A,B) = ^^^{A,B), S{A,B) = {A,-B), T{A,B) = {A,-B). 


Typically, S corresponds to the reflection symmetry in the original equation, when Ginzburg-Landau 
is derived as an amplitude equation, since, say in the Swift-Hohenberg case, u ~ Ae^^ + c.c.. Gomplex 
conjugation symmetry ST is a “normal form artifact. The system is also Hamiltonian, with Ti = 
-|- |Ap — 5|A|^, so that 

A, = dsn, B^ = -dxH. 

One can factor the gauge symmetry and consider the equation in new variables for A,B G C^, the 
Hilbert invariants of the circle action 

a = |Ap ^0, (3 = \B\^ ^0, M = ^{m + \M) = AB G C, with relation 4 q;/3 = MM. 
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The inverse transformation gives A and B up to a relative phase, 


The equations read 


A = B = 


M 

-y/4a/3 


(4.4) 


ax = m 
Px = 'm{a - 1 ) 
nix = 2{/3 + — a) 

Mx = 0. (4.5) 

Since /? is given implicitly through the relation and M is conserved, we end up with a two-dimensional 
system for {a,m), in which we can eliminate (3, most conveniently using the Hamiltonian, which gives 


ax = m 


nix = 2\ 'H — 2a + -a‘ 


(4.6) 


From solutions to this planar ODE, for any fixed we find solutions to the original system by first 
reconstructing (3 from Ti, reconstructing = a/J — from the relation, and using the inverse 


transformation (4.4). Note however that each solution of (4.6) yields two solutions corresponding to 


the different signs of M. 


4.2 Stable manifolds 

Periodic patterns are of the form A = \/l — which in invariant coordinates gives 

a = l-k^, ^ = k^(l-k^), ni = 0, M = -2k(l - k^), B = ^(1 - k^)(l-h 3k^). 


Note that B. \s a. monotone function of within the Eckhaus boundary k"^ G [0; 5)- We can also use 


the expressions for B and M to express (4.6) using the parameter k, 


ax = m 


nix = 2[ B{k) — 2a + j > B{k) = ^(1 “ ^^)(1 + 3A:^) 


(4.7) 


Equilibria are a± = |(1 ± ^1 — \B), or, 


1 


a+ = 1 — a- = -(1 + 3A:^). 

3 

with a+ corresponding to the Eckhaus-stable equilibrium for k'^ < 1/3. We normalize the Hamiltonian 


of (4.7), 


1 

h{a,m-,k) = -m'^ — V{a;k), V{a]k) = / ((1 — A:^)(l + 3A:^) — 4o + 3a^) da = (a — a*)(a — a+)^, 

2 Ji-k^ 
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Figure 4.1: The stable manifold of wave trains, plotted as (a, m, M)(fc, r), with |fc| ^ l/-\/3 (left 3 figures). 
Stable and unstable manifold (right). 


where a* = 2k‘^. The stable manifold of a+ consists of an unbounded branch with 


m 


= —\/2{2'Ha — 2a^ + + /i+), a > a+, 


and a bounded branch, 


m 


= ±\/2(2?^a — 20? + + h+), a* < a < a+. 


A somewhat more compact representation is achieved in coordinates r = , where VF® is given 

through 


a 


= 2A:^ + m = - (1 - 3A:^)), t ^ n = a/2(1 - S/c^). 


(4.8) 


In summary, given a wavenumber |A:| < l/\/3, we can parameterize its stable manifold by r via (4.8). 


Figure 4.1 shows the stable manifold of the wave trains in (a, m, M)-space. 


4.3 Phase shifts 


Phase shifts (p are defined so that for a boundary layer 

lim Im (log(A(x)) — {kx — ip)) = 0, 

X^OQ 

which gives, using the expression for the derivative of A, 

/■“ / M \ 

ip = — lim (Im(logA(x)) — kx) = — Imlog(A(0)) — Aip, Aip = — fed-^ dx. (4.9) 

x^oo Jq \ 2a{x)J 


Using (4.8) gives rr' = m and r' = ^(r —r*)(r+r*) with r* = y^2(1 — Sk"^). We can use r as integration 
variable with r(oo) = —r*, to obtain 


. r^* 2A: , f J2{l-3k^)\ /r(0) 

Aip = — —n-?rdr = arctan ;- + arctan —— 

^ 7,(0) 4A:2 + r2 \ 2k \ 2k 


and, altogether. 


7 V2(l-3fc2) \ ^ /r(0 ) 

ip = arctan | - — - — Imlog(A(0)) + arctan 


2k 


2k 


(4.10) 
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4.4 Computing displacement-strain curves for specific boundary conditions 

General boundary conditions will not be invariant under the gauge symmetry. Consider a general affine 
boundary condition of the form 


B = kiA + k2A + ^, ki,k2, € C] 


(4.11) 


note that this comprises all affine boundary conditions which are not strictly Dirichlet in either com¬ 
ponent. We will discuss Dirichlet boundary conditions at the end of this section. 


The group orbit of the two-dimensional plane defined by (4.11) under the gauge symmetry is 


{iA,By,B = (feiC + k2^ +^i),A = e'-^C C G C,G [0, 27r)}, 
which gives the parameterized representation in invariants 

a = |e|^ M = 2akl^ + k2^+^i). (4.12) 

One now proceeds to find boundary layers by solving (|4.12 ) together with (4.8) and the relation for M, 


a = 2k‘^ + m = “ (1 “ 3A:^)), M = — 2fc(l — fc^), r ^ -\/2(l — S/c^). (4-13) 


Equations (4.12) and (4.13) can be viewed as a polynomial system of 6 real equations in 7 real variables 


Re Im a, m, M,k,T, together with an inequality constraint. We will next show how to solve this 
system in several simple cases. 


Neumann boundary conditions For ki = k 2 = 0, ^ > 0, we can use the second equation in (4.12) 


and the second and third equation in (4.13) to express ^ in terms of r, k, which we can then substitute 


into the first equations of (4.12) and (4.13) to obtain an equation in r, fe, only. After eliminating a 


factor 4A:^ -(- which only gives the trivial solution A = 0, we find 

4/c^ - 8^2 -b (r^ - 2f + Sk^ir"^ - 1) = 0, 


(4.14) 


combined of course with the inequality r ^ -^2(1 — 3k‘^). This equation is a simple quadratic in k = k"^ 
and fj = 

fj -|- (4 — 2\/3 )k — 2 -|- -|- (4 -|- 2\/3)k — 2 — ~ ~ g' 

with solutions given by hyperbolas in a, k, with a double point at \fj,\ = l/\/6 ~ 0.408. 

Solving and substituting back gives the four solution branches which are plotted in the (r, fc)-plane in 


Figure 4.2 



(4.15) 
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Figure 4.2: Plots of solutions (4.15) in the (r, fc)-plane, with reconnection crises. Also shown the restriction 
T ^ r* as a half-ark in the right half-plane. Parameter values are /i = 0.1,0.4,0.41,0.46,0.49,0.7065, left to 
right, top to bottom. 



Figure 4.3: Stable manifold in (a, m, M)-space, together with the Neumann surface and the intersection; pa¬ 
rameters are /r = 0.1,0.4,0.41,0.46,0.49,0.7065, left to right, top to bottom. 


A second crisis happens at ^ = \/2/3 ~ 0.471 when two curves join within the Eckhaus boundary. A 
third and last crisis occurs when /r = l/\/2 ~ 0.707, when the origin becomes a critical point of (4.14) 
and a curve vanishes at the origin. After this crisis, only one curve remains which covers the entire 
range of allowed wavenumbers. 


We can now reconstruct the phase using the original boundary condition A 


^ and (4.10). 


Figure |4.4| shows the resulting displacement-strain curves. For small /x, there are two families of 
boundary layers. One selects wavenumbers, with phase-winding number 1, the other family selects 
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Figure 4.4: Plots of displacement-strain curves in the (</?, fc)-plane corresponding to Figure 4.3 Note that for 
one of the displacement-strain curves has winding number 2. Open ended curves for small ^ correspond 
to curves terminating on r = r,; parameters /r = 0.1,0.4,0.41,0.46,0.49,0.7065, left to right, top to bottom. 


both wavenumber and phase, terminating on r = r*. Along the curve a defect has nucleated from 
the boundary and the defect location diverges to infinity as the end point of the curve is approached. 
At the end point, the phase difference between this boundary layer and the other boundary layer is 
precisely tt, the phase shift along the defect. For increasing /i, the wavenumber-selecting boundary 
layer covers an increasing interval of /c-values, until it reaches the Eckhaus boundary at the extremum, 
at which point it splits into two families of boundary layers, with opposite signs of displacement-strain 
derivative (and opposite stability properties). One of those families eventually reconnects with the 
third family, creating a wavenumber selecting family with winding number 2. 


Pure wavenumber-selection. A special case are boundary conditions of the form 

B = /iA, /i G C, 


(4.16) 


fixed. Since the boundary conditions are invariant under the gauge symmetry, the phase of boundary 
layers is arbitrary once we find a boundary layer. In invariant coordinates, we find the equivalent of 

m — iM = 2A1 = 2^,a, 


a curve rather than the surface we had seen in the previous example. Together with (4.13), we find 6 
equations in 6 variables, which reduce to the system of polynomial equations in r and k 


-2k{l - k^) 

(1-3/c2)^ . 

One readily notices that for Re(/i) = 0, either r = ±r* and k = — Im(/i), the compatible equilibrium, 
or r = 0, Im(//) ± y'T^TpmpI))^. For | Im(^)| < l/\/3, the first ^-solution lies within the Eckhaus 
boundary, for |Im(//)| > l/\/3, the second /c-solution does. 

For small Re(/i), the intersection with the equilibrium splits into a pair of intersections close to ±r*, 
respectively. For Re(/r) > 0, |r| > |r*|, so that only one intersection with the level set of the Hamiltonian 


2Im(;u) (^k^ + = 

2Re(/i) ^2A:^ -|- = 
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Figure 4.5: Plot of selected wavenumber with boundary conditions (4.16) as a function of Im(/i) for fixed 
Re{p) = —.155, 0, .155 from left to right. 





Figure 4.6: Stable manifold in (a, m, M)-space, together with the boundary conditions and intersections. Pure 
wavenumber selection (Re(/r) = —0.1,0.1, Im^ = 0.3, left and middle) and pure phase selection (right). 


yields an intersection with the stable manifold. For Re(^) < 0, both intersections yield boundary layers. 
The dependence of k on Im(/i) is shown in Figure 4.5 The associated pictures in (a, m, M)-space are 


shown in Figure 4.6 


Pure phase selection. Boundary conditions Im(A) = Re(B) = 0 encode the reflection symmetry S. 
In fact, all periodic solntions satisfy this boundary condition for an appropriate phase shift. In order 
to study these boundary conditions more systematically, we write for the group orbit of the boundary 
condition A = B = which gives 

m = 0, M = —2i(^ — ^^). 


Solving the remaining equations (4.13), we find 4 solutions corresponding to the exact periodic pattern, 
r = ±r*, ^ = ±(1 + ik)V^ — and two solutions r = 0,^ = (i(l — k‘^)/2 — k)jy/2. The first four 
solutions correspond to only two actual “boundary layers” (r = ±r* yield the same boundary layer). 
They are in fact the exact periodic solution and the solution with a phase shift of tt. The latter two 
solutions correspond to an appropriate shift of the defect and yield a phase shift of 7r/4. 


A completely equivalent analysis also yields the boundary layers for the third reflection symmetry 
Re A = 0, ImR = 0, with the same wavenumber selection, and phases shifted by 7r/4. 
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Dirichlet boundary conditions. For Dirichlet boundary conditions, A = 
form gives 


{(^, B); B = A = e'V, ? G C, G [0, 27r)}, 


/i > 0, the group invariant 


which gives the equivalent of (4.12), 


a = ^ 


M. = 


(4.17) 


Again, we can use the second equation in (4.17) and the second and third equation in (4.13) to express 
^ in terms of T,k. The first equations of (4.12) and (4.13) give 


1 

2 


+ 2/c^ = /i 


2 . 


(4.18) 


the solution curves, together with the restriction r ^ r* are shown in Figure 4.7 The ellipse intersects 
the boundary r = T^{k) for /i > 


As \x is increased, we see displacement-strain curves changing from wavenumber selection (winding 
number 1) to phase-selection. The second curve of boundary layers terminates at defects and, as [i 
increases further, shrinks and vanishes on the defect associated with k = 0. 


Note that for fj. = 1.1, ^{k) is not monotone. In fact, <p{k) develops 3 critical points for /x < 1 in a 
pitchfork bifurcation. This phenomenon cannot occur for well-posed equations in the phase-diffusion 
systems, where an eigenvalue or resonance pole disappears at infinity at ^p' = 0. 






Figure 4.7: Plots of solutions (4.18) in the (t, fc)-plane, with reconnection crises. Also shown the restriction 
T ^ T* as a half-ark in the right half-plane. Parameter values are /x = 0.2,0.7, 0.93,1.1. 


Remark 4.1 All boundary layers constructed in this section are in fact transverse in the sense of 
Definition 2-4 1 provided that the intersection in the reduced orbit space is transverse. An example of a 
non-transverse boundary layer is the half-defect A^ = itanh(x/\/2)g, x > 0, with boundary conditions, 
Ax = Inspecting the formulas for the k-dependent defects (4-^) o,nd differentiating at the origin 

X = 0, we find tangent vectors [A, B) = (0, —1) from the phase, {A, B) = (i, 0) from x-translation, and 
{A, B) = (2,0) from the k-derivative, so that the real span of (0, i) is not contained in -\-TB. The 

reason is of course that these boundary conditions mimic infinitesimal translation of the defect at the 
origin. 
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Figure 4.8: Stable manifold in {a, m, M)-space, together with the Dirichlet plane and the intersection. Parameter 
values are fi = 0.2,0.7,0.8,0.82,0.93,1.1. 







4.8 


We first see 


Figure 4.9: Plots of displacement-strain curves in the {ip, fc)-plane, cporresponding to Figure 
the separation crisis at the Eckhaus boundary, followed by the vanishing of two defect limits. Eventually, p is 
not monotone. 


4.5 Existence — linear and afRne boundary conditions 

Linear boundary conditions translate into a cone in a, m, M-space, by scaling invariance. On the other 
hand, the expansion of VF® at a = m = M = 0is quadratic, so that any cone is contained in the interior 
of the body of the fish, which guarantees intersections between VF® and the cone, hence existence of 
boundary layers. 


22 
































































Proposition 4.2 For all affine boundary conditions 

TiA + T 2 B = /i, , 

with Tj : —)■ real linear, /U G RgTi + RgT 2 = there exists a boundary layer. 

Proof. The proof exploits the topology of (a, m, M)-space. We can distinguish three different cases 
depending on the rank of T 2 : 

(1) rkr2 = 2; 

(2) rkTs = 1; 

(3) rkr2 = 0. 

In the following, we will treat various scenarios on a case-by-case basis. 

Case (3). The map Ti is invertible, so that the boundary conditions are equivalent to Dirichlet 
boundary conditions, and we already established existence in the previous section. 

Case (1). We can invert T 2 and obtain a boundary condition of the form 

B = TA + V, T = Tf^Ti, V = 

or B = kiA + K 2 A + V. In case ki = K 2 = 0, Ti = 0, we recover the Neumann boundary conditions 
studied in the previous section and we can conclude existence. We consider the boundary manifold in 
the space of invariants a, m, M, thinking of m as the vertical coordinate. We set A = f = rF^ , which 
gives 

a = r^, M. = r"^ -|- K2e‘^^'^^ + ruF'^. 

We view this as a family in r of closed curves, each parameterized by and compare the location 
of the curve with the position of the stable manifold. For large r, the stable manifold is located near 

m ~ —a^, so that the curves M. ^ r"^ (/ii -|- K 2 e^'’^) lie above the stable manifold. For small r, suppose 

first that v > 0. The curve M ~ ruF^ then is located either inside the compact part of the stable 
manifold, or it winds around or intersects it. In either case, homotoping between small and large r, 
we see that there necessarily is an intersection. In the case /x = 0, linear boundary conditions, the 
small-r asymptotics are equal to the large-r case and one readily sees that the curve is located inside 
the compact part of the stable manifold, which has a quadratic tangency with the plane a = 0. Again, 
a homotopy to large r gives intersections between our family of curves and the stable manifold. 

Case (2). Without loss of generally, possibly rotating with the gauge symmetry, the boundary con¬ 
ditions can be written in the form 

Re 5 = Re A -|- K 2 Im A -|- , Im A = cji Re A -|- ;U 2 , 

with real parameters ^i,ii 2 ,ki,H 2 ,(Ti. Parameterizing ImB = y, Re A = x, we find 

a = + (cJix -I- ^2)^, M = {x + i(cJix -|- ^2)) + K2crix + Hi + K2li2 + iy) • 
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We view this surface as a family in x of straight lines parameterized by y. We follow the family of 
lines from x = +oo to x = —oo, and note that near the endpoints of the homotopy, the lines are close, 
with tangent close to (i — (Ti)x^. Since M. = 0(a) when y = 0(x), the line is located above the stable 
manifold for large |x| (and large a), with non-real tangent. As x is increased from —oo to -|-oo, the 
tangent vector rotates by tt. As a consequence, the family of lines necessarily intersects the stable 
manifold. ■ 


In the case where the equation comes from a variational problem, 

— J ~ dx-|- (/r|Ap -|- (i/A^ -|- c.c.) -|- {pA + c.c.)) , /x G M, p G C, 

one obtains boundary conditions 


B — pA. vA Pf P ^ P ^ 

which form a Lagrangian (affine) subspace of C^, with respect to the symplectic form 

^{{AiiBi), {A2, B2)) = Re(Ai.B 2 — A2B1). 


(4.19) 


Since in the interior of the domain, | A| = 1 and the associated solutions A = e'*^ are the only minimizers, 
one expects boundary layers with finite energy, thus connecting to /c = 0. While such a result can 
presumably be established using variational techniques, possibly for more general systems, we restrict 


here to the Ginzburg-Landau equation and the geometric viewpoint from Proposition 4.2 


Proposition 4.3 For all Lagrangian boundary conditions (4-19), there exists a boundary layer with 
k = 0. 


Proof. We describe the boundary surface in the invariant space setting A = ^ = re'^, so that 

a = r"^, M = {p + i>e^''^)r^ -|- pF'^r. 

For p = 0, these boundary conditions describe a cone with elliptical cross-section, which necessarily 
intersect the compact part of the stable manifold in the plane M = 0. For p A 0; we can view the 
boundary manifold as a family in r of curves parameterized by ip. One can choose two intersection of 
these curves with the real axis (maximum and minimum) depending on r in a continuous fashion (the 
curves are Limagons). For r large, this intersection point lies above W®. For r small, the intersection 
points lie either above and below, or inside, or on the compact part of the stable manifold. For large r, 
the intersection points lie above the stable manifold. As a consequence, there is an intersection, which 
concludes the proof. ■ 

Remark 4.4 (Non-existence and isolas) Boundary layers may not exist even with well-posed bound¬ 
ary conditions. A simple and instructive example is the cubic-quintic Ginzburg-Landau equation 

At = - A + 7 A|A |2 - A|A|^ 

Periodic orbits re'^^ exist for k‘^ = — 1 -t- provided 7 ^ 2 . The steady-state equation preserves 

the Hamiltonian H = \Ax\‘^ — |Ap -|- ^jA]^ — ||A|®. For 7 ~ 2, all periodic orbits are close to |A| = 1 
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and have H ~ —1/3. Choosing “asymptotic boundary conditions”, B = W^{A = t),B = 0), we 
conclude that all boundary layers would have Hamiltonian H = t), which is incompatible with the 
Hamiltonian of periodic orbits. The unstable manifold W”" can in fact be expressed as a smooth graph 
B = so that the equation with these nonlinear boundary conditions is in fact 

locally well-posed. Alternatively, imposing Dirichlet boundary conditions A = 0 gives H ^ 0 on the 
boundary, which also precludes intersections with the stable manifold of periodic patterns, where H ~ 
-1/3. 

Considering inhomogeneous Dirichlet boundary conditions A = p. with increasing pL, one finds a tan- 
gency between boundary conditions and stable manifold (much like the tangency in the cubic case at 
p. = 0, just at some finite a > 0, k = 0), which then leads to a small closed curve of intersections and 
a small isola in the (</?, k)-plane. 

4.6 Scaling boundary conditions 

In general, Ginzburg-Landau spatial dynamics are recovered on a 4-dimensional center manifold for 
spatial dynamics [ 13 ]. Boundary conditions will typically intersect the center-stable manifold along 
a two-dimensional manifold,which can be collapsed along the stable hbration to obtain effective 2- 
dimensional boundary manifolds in the center manifold. In the center manifold, the dynamics are, up 
to changes of coordinates and higher-order terms. 

Ax = ikA -\- B, Bx = ikB — pA 

where p is the bifurcation parameter. The usual scaling A = y/pA, B = pB, k = y/pk, y = y/px, 
eliminates the parameter p and gives, up to the fast rotation, the Ginzburg-Landau equation. Under 
the scaling, most linear boundary conditions will reduce to homogeneous Dirichlet boundary conditions 
at leading order since derivatives are small in the expansion. 


5 Computation of boundary layers in the Swift-Hohenberg equation 


Computing boundary layers is related to the computation of heteroclinic orbits. Numerically one needs 
to approximate the infinite domain [0, oo) by a hnite domain [0, L] and discretize in space. The first 
step, truncation, involves choosing appropriate boundary conditions at x = L. The correct boundary 
condition, u{L) G is usually approximated by its linearization. A difficulty in our case is the 
presence of neutral directions, phase and wavenumber, and the fact that the asymptotic state is not 
explicitly known. 


In the following we present the general strategy of our approach, Section 5.1 and some results for a 


specihc set of boundary conditions. Section 5.2 


5.1 Decomposing periodic and heteroclinic orbit 

Let x(x) be a smooth function with x(x) = 0 on x < I’-, x{x) = 1 on x > £+, for some 0 < ^_ < £+ < L. 
Also, denote by = —{dxx +1)^ + T the linear part of the Swift-Hohenberg equation. Lastly, we write 
X]u := C^{xu) - x^u^- 
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With these preparations, we substitute the ansatz u(x) = x{x)ustikx — (p]k) + w{x) into the Swift- 
Hohenberg equation, subtract x{Lust — Ugt)^ obtain 

^o{w, ip, k) = - {{xust + wf - ixustf) + [£^, - (x - X^)^^st = 0> 

{w,Wx,WxxiWxxx)iS^) G B. (^'l) 


Since x vanishes near x = 0, w needs to satisfy the same boundary conditions that we imposed on u. 
The hrst two terms vanish for tc = 0, reflecting the fact that Ugt solves the Swift-Hohenberg equation. 
The last two terms, which we refer to as residuals, are compactly supported on [.£_,£+]. 

We are interested in solutions w that are exponentially localized. Define therefore = {tt; u{x)e^^ G 
L^}, with induced norm, and, analogously = {u]d^u G ^ k}. It turns out that Tq is not 

differentiable in with respect to /c as a map: the derivative with respect to k of the term 3(xust)^u; does 
in general not belong to due to the linear multiplier Bx^UstxdxUst- A simple rescaling can remedy 
this difficulty as we shall see below. Consider the scaled differential operator = ~{k'^9yy + 1)^ + 
with appropriately scaled boundary conditions, and define 

To{w, if, k) = - {ixust + wf - ixustf) + x]^^st - (x - = 0) 

{w, kWy, k^Wyy, k^Wyyy){^) G B, (5.2) 


where u^t = Ust(y “ fik) is the 27r-periodic, scaled family of stripes, x = xiu) can be left unsealed. 
One readily finds that solutions w,k, ip to (5.2) are in one-to-one correspondence with solutions w,k,ip 
to (5.1) via the rescaling w{x) = w{kx) and the appropriate modification of x- 


Lemma 5.1 (Well-Posedness) The map J-q : n S x —)• is well-defined and smooth near 

a given boundary layer w = Ubi — X^^st? for rj > 0, suffieiently small. Moreover, the linearization at a 
boundary layer is Fredholm with index 1. If in addition the boundary layer is transverse, according to 
Definition 2-4, then J-'q is onto and the k- and ip-eomponents of the one-dimensional kernel do not both 
vanish. 


Proof. The map is well-defined since residuals are compactly supported, thus belong to exponentially 
weighted spaces. The linearization with respect to w is Fredholm with index -1 as a simple counting 
of Floquet exponents shows. Adding two parameters increases the Fredholm index to 1. The map 
and derivatives with respect to w depend smoothly on k and ip as coefficients ttgt depend smoothly 
on k and ip in L°° and the dependence of the linear operator on k is smooth. The last claim is an 
immediate consequence of the assumption of the intersection between boundary manifold and center- 
stable manifold in Definition 12.41 ■ 

The Lemma can be viewed as a substitute for the smooth dependence of the strong-stable foliation of 
the periodic on phase and wavenumber, parameterizing the base points of the fibration. As a direct 
consequence of the lemma, we obtain the existence of a one-parameter family of boundary layers, 
described by a smooth curve in the [k, (^)-plane. Using Lyapunov-Schmidt reduction. 

Truncating to a finite domain would require w{L) to he in the (unknown) strong stable fiber of Ust 
(translated by ttst)- Since w is exponentially small, this fiber is well approximated by its tangent space, 
so that a good asymptotic boundary condition would be w{L) G E^^{kL — ip\k), where E^^{y) denotes 
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the strong stable Floquet subspace of the linearization at Ust{y)- Again, we find the computation of 
this Floquet subspace somewhat intricate. Alternatively, one can impose other boundary conditions 
at X = L, as long as the boundary subspace is transverse to FI™. In practice, we choose Dirichlet 
boundary conditions and add a phase condition, to achieve the correct number of boundary conditions. 
As a phase condition, we choose 

I'K-\-27V 

Tph{w,ip)= u's^iy - (p)w{y)dy, (5.3) 

Jk 


which ensures that w cannot encode a phase shift in the periodic pattern (which is already encoded 
explicitly in the parameter cp). 

We obtain the following problem in a bounded domain: 


Tq{w, if, k) = - {{xust + wf - ixustf) + x]nst - (x - X^)^^st = 0: 

{w, kWy, k^Wyy, k^Wyyy){0) G B 
{w,Wyy){L) = 0 

rK-\-27T 

u'^^iy - ip)w{y) = 0. (5.4) 

Jk 


One readily verifies that this problem is Fredholm index 1 and we therefore expect to find solution 
curves in the (F, y5)-plane. We expect exponential convergence of the solutions to (5.4) to solutions of 
(5.2), exponentially as L —)• oo. 


5.2 Wavenumber selection in Swift-Hohenberg — numerical results 


We use second order finite differences to discretize the fourth-order differential operator, using a mesh 
of width /i on a domain of length L. As part of the problem, we compute Ust{y\k) using Neumann 
boundary conditions on a domain of size 27r. We then construct Ust{y — 'P',k) using cubic splines and 
periodic extrapolation. We then solve (5.4) using arc length continuation within Matlab. We found 
the expected exponential convergence in the domain length, which breaks down near the Eckhaus 
boundary. We verified that results are independent of the cut-off function. 


In order to illustrate the phenomena discussed here, we computed the displacement-strain curves for 
free boundary conditions, 

'^xx An — 0, '^xxx A Ux — 0. (5.5) 


Free boundary conditions arise naturally when considering the operator generated by the bilinear form 
f^{dxxU A u){dxxV A n)dx on x H^. Figure 5.1 shows results of these computations for various 


values of y,. Numerical parameters were dx = 0.04, L = 14 * vr; for these parameters, we found errors 
in the displacement-strain relations below 10“^. 


Not surprisingly, we hnd phase selection for small //, reflecting the fact that the energy density along 
periodic patterns is almost constant in x. For moderate sizes of y, we still find wavenumber selection, 

^Such an approach is well understood in the context of heteroclinic or homoclinic orbits to hyperbolic equilibria. One 
finds that the convergence of solutions is with correct linear Floquet boundary conditions and only e“’’^ with 

general transverse boundary conditions, where rj measures the spectral gap separating Floquet exponents on either side 
of the imaginary axis; see m §6] and references therein. 
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/i =0.01 /i =0.1 /i =0.3 /i =0.8 ^=1.5 /i =1.9 ^=2.02 /i =2.05 ^=2.5 /z =4 


Figure 5.1: Displacement-strain relation in the (fc, :^)-plane for in the Swift-Hohenberg equation with boundary 
condition (5.5), for various values of /i, fc S (0.6,1.1). The vertical lines indicate Eckhaus boundaries. Note the 
separation crisis at the Eckhaus boundary near p. = 2.02 followed by the reconnection with a different branch 
(near /r = 2.15). 





/i =2.5 

Eigure 5.2: Displacement-strain relations as in Eigure 
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k = 0.91 
V? = 3.00 


k = 0.93 
If = 3.66 


k = 0.9276 
^ = 11.31 
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5.1 


for fi = 2.5,4, extended for large phases (left). Plots 


of associated profiles for /r = 4 associated with points marked on the displacement-strain curve for /r = 4; 
increasing top to bottom, left to right. One clearly notices the separation of an interface from the boundary 
along the snaking curve (right). 


with a characteristic double-well shape for k > 1 that persists throughout the entire range of Note 
also that the displacement-strain curves are symmetric by i—)> (/?-|-vr but break the (p i ——p symmetry. 
For ~ 2.02, the displacement-strain curve touches the Eckhaus boundary and separates, the typical 
transition between wavenumber and phase selection. Shortly after, /r ~ 2.15, the displacement-strain 
relation reconnects with a different branch, which is not periodic in p, but rather “snakes” as (/? —)> oo. 
The profiles indicate that this snaking behavior is mediated by the emergence of a stationary interface 
between the stable state u = — 1 and the periodic patterns. This stationary interface interacts 

with the boundary condition in an oscillatory, exponentially weak fashion, yielding steady states with 
increasing separation of the interface from the boundary and exponentially small, oscillatory oscillations 
in the wavenumber. 


6 Discussion 

Summarizing, we have introduced a systematic description of boundary layers in terms of curves that 
can be interpreted as displacement-strain relations. Generically, in terms of the boundary conditions, 
these are smooth curves which may terminate either near instability boundaries, or when the boundary 
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layer itself disappears. We showed how such curves can be used to describe systematically the set of 
equilibria in bounded domains, in particular as the domain size is varied. More specifically, we defined 
an operation on pairs of displacement-strain relations corresponding to left and right boundary layers 
in large but finite domains, that consists of “vertical subtraction” and subsequent “quantization”. 
The resulting families of intersection points are in one-to-one correspondence with equilibria in large 
bounded domains. 

Most of the paper is devoted to actually constructing or computing such displacement-strain curves. 
We saw that the common phase-diffusion approximation only gives very limited classes of examples! On 
the other hand, the Ginzburg-Landau equation, equipped with a variety of inhomogeneous Dirichlet, 
Neumann, or mixed boundary conditions gives a plethora of interesting phenomena, including phase- 
and wavenumber selection, non-monotone displacement-strain, pinch-off of defects leading to end points 
of displacement-strain curves, absence of boundary layers and isolas, and displacement-strain curves of 
higher winding number. We also outline how to systematically explore such displacement-strain curves 
using numerical continuation and illustrate some preliminary results in the Swift-Hohenberg equation. 

Our point of view raises numerous questions. First, we noted throughout that changes in the type 
of displacement-strain curve are accompanied by reconnection or separation events at the Eckhaus 
boundary. A local analysis of the stable manifold there should reveal a universal description of such 
events. One may also wonder if a variational structure in the boundary condition may lead to additional 
restrictions on boundary layers and bifurcations.We showed existence in the Ginzburg-Landau equation, 
but we suspect that, more generally, variational problems where periodic structures minimize the energy 
always accommodate boundary layers. 

Much of our analysis is reminiscent of the analysis of snaking diagrams in, for instance, the weakly 
subcritical Swift-Hohenberg [T]. There, the trivial stable state coexists with stable periodic patterns 
and interfaces between trivial and periodic states provide “effective boundary conditions”. Due to 
translation invariance, such boundary conditions are always pure wavenumber selecting, k{(p) = const, 
so that interesting questions are concerned with changes of the selected wavenumber as parameters 
/i are varied. Fold points of /c(/r) can be interpreted as unpinning transition. The analysis in [T] is 
concerned with phase matching of such “effective boundary layers”, which yields spatially localized 
structures. Due to the singularities of displacement-strain curves, equilibria are not locally unique, 
generating various bifurcations (ladders, etc.). It would clearly be very interesting to analyze phase 
matching as described here near singularities of displacement- strain curves, and possibly elucidate the 
connection with snaking diagrams further. 

In higher dimensions, rolls allow for rotation as an additional degree of freedom. Setting up boundary 
layer problems in x > 0, and y G 'RjLX periodic, we do encounter however a translational symmetry 
in y which implies wavenumber selection due to arbitrary phases, for rolls that are not parallel to the 
boundary. The case of Neumann boundary conditions has been studied extensively, since boundary 
layers can then be interpreted as symmetric grain boundaries; see [EHH and references therein. 

Recognizing the analogy with defects such as grain boundaries, one can ask if there is a similar de¬ 
scription for defects and their interactions. This appears to be an interesting avenue of research. Some 
results in two space-dimensions [7] indicate that properties such as phase selection can be attributed 
to inhomogeneities in higher-dimensional striped phases. 

Generalizations of our setup also include discrete atomic lattices, possibly multi-atom lattices, which 
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possess natural periodic equilibria. In a different direction, boundary conditions (and defects) for 
oscillatory media, where striped patterns propagate with non-zero group velocities, have been classified 
in [12]. It would be interesting to find a more unified description of all these interesting phenomena, 
and we hope that the present work is a contribution in this direction. 


Acknowledgments. The authors acknowledge discussions with David Lloyd on the implementation 


of the computational strategy outlined in Section 5.1 The decomposition described there was first 


implemented in a joint forthcoming work on grain boundaries [9| and the computational results in 
Section 5.2 mimic the approach in |^. 
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